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Abstract 

We present a semiclassical and mostly analytical model of elas- 
tic neutron breakup reactions for exotic projectiles. Both nuclear 
and Coulomb induced reactions are considered and the potentials are 
treated to all orders in the interactions. Furthermore we introduce 
a technique which allows the use of the full Coulomb potential, thus 
including all multipoles besides the dipole. Results for deeply bound 
states as well as for halo initial states are presented and it is shown 
that experiments on heavy targets would be well suited to study exotic 
nuclei with tightly bound valence nucleons. 

Pacs 21.10.Jx, 24.10.-i, 25.60.Gc, 27.30.+t 



1 Introduction 

The usefulness of breakup reactions in addressing the problem of exotic nuclei 
structure is widely recognized since many years. In particular, the study of 
spectroscopy and single-particle strength of halo nuclei allows a direct com- 
parison between theory and experiment, but in order to do so reliable reaction 
models are needed to obtain accurate cross-section calculations. Therefore a 
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good understanding of the reaction mechanism is an issue of capital impor- 
tance. 

This work regards one-neutron breakup from exotic nuclei, with both 
halo and tightly bound valence particles. We concentrate on reactions on 
heavy targets at intermediate beam energies. The projectile ground state is 
regarded as a single-particle state in the two-body system made of a valence 
neutron plus a heavy fragment. Most of the experiments are inclusive and 
only the heavy fragment is detected, but sometimes exclusive experiments 
with the valence neutron detected in coincidence PP-|3] with the core have 
been performed. Since we consider reactions involving heavy targets, the 
process is Coulomb dominated for halo projectiles. However we will show 
that for normally bound initial states Coulomb breakup becomes negligible 
and breakup reactions can be used as a tool to obtain structure information 
provided nuclear breakup is properly included. 

There have been several approaches to this problem, some of which 

consist in numerical solutions of the Schrodinger equation [22J- 27J. These 
approaches have reached a large degree of accuracy but the physics inter- 
pretation of the results partially relies on the comparison with results from 
more analytical models. Furthermore due to the numerical complexity the 
scattering of two heavy nuclei has not been treated so far. Thus our focus is 
not so much on the numerical aspects but rather on the understanding of the 
main features of the reaction mechanism in a simple, yet accurate, picture. 
In this spirit, earlier works |2*H1 12*9] have shown the importance of treating 
the nuclear and Coulomb interactions in a consistent way and to all orders. 

Within the semiclassical formalism, we now attempt to examine two as- 
pects of the reaction process: the first one is the peripherality of the reaction, 
whereas the second concerns the accuracy of the dipole approximation for the 
Coulomb interaction. They are both important to extend breakup studies 
to heavy reaction partners. The concept of the reaction being peripheral is 
important as it allows the use of spectator-like models for the core target 
interaction leading to a semiclassical relative motion trajectory. It underlies 
a large part of existing approaches, and has been recently rediscussed by 
Capel et al.j2Zj. The fact that the heavy fragment is detected, i.e. it sur- 
vives the collision suggests that processes where it collides closely with the 
target would not contribute significantly to the measured cross-section. The 
main contribution is expected to come from interactions between the target 
and the tail of the valence neutron single-particle wave function. In halo 
nuclei this tail goes much beyond the range of the potential. One of the aims 
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of this work is to compare the cross-sections obtained by using the whole 
single-particle wave function with that given just by its tail. This turns the 
formalism of [2E1 I2H] less analytical and adds some numerical complications, 
but some light can be shed on the problem of which part of the projectile 
wave function the reaction is sensitive to. This is an interesting point being 
at the center of a scientific debate since recently the reliability of spectro- 
scopic factors obtained from peripheral reactions has been questioned |3Uj . 
For deeply bound states peripherality is best obtained on heavy targets. But 
as we will show in the following in these cases Coulomb and nuclear breakup 
can be of the same order depending on the nucleon separation energy and 
they need to be treated consistently up to separation energies around 8 MeV. 
However for larger separation energies and angular momentum states t > 1 
Coulomb breakup becomes negligible. Therefore for reactions involving very 
neutron rich projectiles we recommend the use of heavy targets. 

Regarding our second aim, the formalism developed in [2E1 I2H] included 
Coulomb potential only in its dipole form, although all orders in pertur- 
bation theory were considered. It was shown that full description of the 
time-evolution could be obtained analytically. On the other hand, other pa- 
pers jTH] have stressed the importance of going beyond the dipole term, in 
particular when one deals with proton haloes. A series of works has also 
focused on the possible effects of dipole-quadrupole interference (20] • The 
second goal of this paper is therefore to calculate the breakup due to the 
full Coulomb potential and to compare the results to those from the sim- 
ple dipole approximation. This is done in preparation for a future study of 
proton breakup. 

The new aspects of the present approach with respect to our previous 
works [2H1 US! are therefore: i) the use of the full numerical Woods-Saxon 
wave function for the initial neutron state: ii) an all-multipoles-all-order 
treatment of the Coulomb potential; iii) the application to reactions in which 
both projectile and target are heavy nuclei. 

2 Theoretical framework 

In previous papers [2E1 121] some of us introduced a full description of the 
treatment of the scattering equation for a projectile which decays by single 
neutron breakup due to its interaction with the target. There it was shown 
that within the semiclassical approach for the projectile-target relative mo- 
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Figure 1: Reference frame and coordinates used in this work. 

tion, the amplitude for a transition from a nucleon bound state ipi in the 
projectile to a final continuum state ipf is given by 

1 f°° 

9fi=~ dt< ^ f (f,t)\V(f,R(t)M(r,t) >, (1) 

^ J —oo 

where V is the neutron-target interaction, responsible for the neutron tran- 
sition to the continuum. R(t) is the core-target relative coordinate, that is 
approximated by a straight line trajectory with constant velocity v , R(t) = 
xR± + zvt, and f is the neutron-core coordinate. As usual in a semiclassical 
formalism, the concept of trajectory brings on the time dependence. Explic- 
itly, the initial state wave function depends on time as if>i(r,t) = e- i£ot/h (l)i(r), 
where Eq is the binding energy of the valence neutron. The final state wave 
function, will be specified later. For light targets the recoil effect due to 
the projectile-target Coulomb potential is rather small and the interaction 
responsible for the reaction is mainly the neutron-target nuclear potential. 
In the case of heavy targets the dominant reaction mechanism is Coulomb 
breakup. The Coulomb force does not act directly on the neutron but it 
affects it only indirectly by causing the recoil of the core. In Ref. [2B] it was 
shown that the combined effect of the nuclear and Coulomb interactions to 
all orders can be taken into account by using the potential V — V nt + V e s 
sum of the neutron-target optical potential and the effective Coulomb dipole 
potential. In Ref. j^EL as in this paper we neglect the neutron-core final state 
interaction under the hypothesis that the observables calculated do not de- 
pend significantly on it (HI]- If for the neutron final continuum wave function 



4 



we take a distorted wave of the eikonal-type 

ipf{r,t) = exp jifc • r + ixeik(r,t) - ( 2 ) 

where k = (k x ,k y ,k z ) is a real vector and the eikonal phase shift is simply 

1 r°° 

Xeik(r,*) = -y \/(f,i?(tW. (3) 
The coordinate system used in this work is presented in Fig^ where 

b = R ± -r ± = R ± - (4) 

Then the amplitude in the projectile reference frame becomes 

9fi(k, b) = ^jd z fj dte-^+^-^^V (V, R{tf) ^ (f) (5) 

where (pi is the time independent part of the neutron initial wave function 
and i = stands for the angular momentum quantum numbers, u = 

(ef — Eq) /h and Eq is the neutron initial bound state energy while £f is the 
final neutron-core continuum energy. This formalism is closely related to 
[3~2] . Eq. (j5J) is appropriate to calculate the coincidence cross section A p — > 
(A p — 1) + n. Finally the differential probability to the neutron momentum 
can be written as 

and we have averaged over the neutron initial magnetic substates. 



2.1 Coulomb phase 

The dipole approximation to the Coulomb potential is 

Vcouir, R) = -J^— - ^ - V fr££ + 0(r\ (7) 
\R - (3 x r\ R R 

where Vq = Z c Z t e 2 and (3\ = m c /m p . The centre of mass of the projectile is 
assumed to follow a straight line, R(t) = xR± + zvt where v is the relative 
motion velocity at the distance of closest approach. 
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As stated above, one of the aims of this paper is to derive analytical 
formulae for the breakup amplitude when the above approximation is not 
performed. As done in j2H], this implies the calculation of the integral 



Xeffib, r,k) = ^J dte iuJt V C ou(r, t), (8) 

with uj = (ef — ei)/h. Such integral is solvable in the dipole approximation. In 
order to obtain a similar result for the whole multipole expansion, a screening 
term is added and subtracted to the potential, thus it can be written as 

V(r, R) = V ah (f, R) + V lo (r, R), (9) 

where 

e --y\R-Pir\ e -yR 
V sh (f, R) = V - V —-, (10) 

\R — pir\ 



1 _ e -j\R-Pir\ 1 _ e -jR 

V lo (f,R) = -V + V . (11) 

\R-/3if\ R 

The term V s h contains the singularity at R = but decays quickly with the 
impact parameter. On the other hand, Vi Q , well-behaved in the origin, ac- 
counts for the long-range character of the Coulomb potential. When inserted 
in eq. (JBJ), these two terms can be treated in different ways if the parameter 
7 is big enough. In this done in |2H] with the nuclear potential, V s h 

can be considered in the sudden approximation, yielding a phase 

Xsudd{b, r) = ~ J dtV sh {f, t), (12) 

whereas V\ needs to keep the whole time evolution description, but, being 
weak, it can be approximated to first order 

X per t(b, r,k) = l [ dte^V lo (f, t). (13) 



Therefore the Coulomb phase becomes a sum of two terms 

Xeff(b, r, k) = Xsudd(b } r) + x P ert(b, r, k), (14) 
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both of them depending upon the screening parameter 7, as 

Xs Udd (b, R) = ^ (#0(76) - MiR±)) , (15) 

X per t(b, f, k) = -JLe^ v (K {bu/v) - Koiy/ibu/vy + ^bf) 

- ^ (k (R ± u/v) - K (y/(R ± u/v)* + (iR±) 2 )) . (16) 

In order for this approximation to be valid, the screening term 7 needs to 
be sufficiently large as to ensure that the range of V s h remains short enough, 
and that V\ does not become too large. This is next achieved by taking just 
7 = 00, in which case 

Xsudd = (17) 

Xpert = -j£ (e^^Koiub/v) - K (uRjv)) (18) 
If eq. (fTSj) is expanded up to first order in r, we obtain 

Xpert ~ ^ ( K (u;R ± /v)— i 3 1 + Kl (u;R ± /v)^^-i3 1 j . (10) 

consistent with the result of 

2.2 Sudden limit and all-order treatment 

Eq. (|18|) gives rise to a probability amplitude which contains only first order 
in the interaction potential Vq. Aiming for an all-order formalism, in j2H] it 
was shown that a possible way to achieve this is to use sudden approximation, 
subtract the first order term, which diverges for large impact parameter, and 
then to add a first order term calculated in time- dependent perturbation 
theory, i.e. eq. (JTHJ). 

The sudden limit (u — > 0) must be therefore taken in the above expression 
for Xeff, yielding 

*5? = ^*±. (») 
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The Coulomb breakup probability amplitude is therefore 

,Cou _ f j^Ak-r, /-* ( JxT/f 1 2V ° - 6 



0<"» = y dre^^f) ^/ - 1 - log — + Xper4 j (21) 
Moreover, the nuclear part of the breakup amplitude is given by 

g nuc = j dfe ik-T ( e i X nt(b) _ ^ (22) 



With these ingredients, the expression for the differential cross-section is just 
da 1 mk 



de An 2 h 2 



j dbb\S ct (b)\ 2 \g nuc + g Cou \ 2 } (23) 



where S ct (b) is the core-target S-matrix, depending upon its appropriate im- 
pact parameter b. Then we can apply this expressions of all-multipoles and 
whole wave function in two formalisms: 

a) - [2H1 Nuclear in sudden approximation and Coulomb to first order in 

Xeff- 

b) - [211 Nuclear and Coulomb in sudden approximation, first order sud- 
den Coulomb subtracted and first order perturbation added. 



3 Applications 

3.1 Two-body interactions 

In the calculations below, the ground state of n Be is assumed to be a purely 
single particle s-state of the two-body 10 Be-n system. Possible deformations 
of the 10 Be core are not considered, thus its intrinsic spin is just zero in our 
approach. The ground state wave function is calculated in a Wood-Saxon 
potential with fixed geometry, radius 1.25A C fm and diffuseness 0.6 fm. 

Our formalism only considers the interaction between the core and the 
target inside the core probability, given by the modulus square of the S- 
matrix S ci . This S-matrix is calculated by making use of a potential V ct 
obtained by double-folding the matter densities of 10 Be and 208 Pb with an 
effective nucleon-nucleon interaction |33j. The densities where assumed to 
have gaussian shapes, with parameters chosen to fit the experimentally de- 
termined root mean square radii, which are 2.30 fm for 10 Be [34 and 5.6 fm 
for 208 Pb jSHl- 
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Figure 2: Calculations of differential cross-section for Coulomb breakup of 
11 Be on 208 Pb with respect to 10 Be-n relative energy in the centre of mass of 
the projectile. Beam energy is 70 A. MeV. Full line: full-range bound state 
wave function. Dashed line: asymptotics. 

Finally, when nuclear breakup is addressed below, a potential interaction 
V nt between the valence particle and the target is required. The optical 
potential of jSHl HH] is used in our calculations. 

4 Spectrum of elastic breakup of n Be on 208 Pb 
4.1 Validity of asymptotics 

This formalism allows us therefore to directly compare calculations involving 
the full-range bound state wave function, calculated as shown in Section ETT1 
to those obtained with the asymptotic form of such wave function, as done 
previously [2B]- Fig. El shows calculations of 10 Be-n relative energy spectrum 
after Coulomb breakup of n Be against a 208 Pb target at an incident energy 
of 70 A. MeV [3]. The asymptotic approximation is indeed very close to 
the full-range wave function, and fully within experimental error bars, as 
will be shown below. Besides, the asymptotics yields a slightly larger value 
of the differential cross-section in particular for the higher values of relative 
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Figure 3: Calculated impact parameter profile in Coulomb breakup of Be on 
208 Pb at 70 A.MeV. Full line: full-range bound state wave function. Dashed 
line: asymptotics. 

energy, as expected. The difference between the two approaches can be easily 
understood from a simple Hulthen-like argument. It must be pointed out, 
however, that this check can only be carried out for s-waves, because, in this 
formalism, the probability amplitude of eq. (}2*2*|) diverges upon integration 
over k, when the asymptotic expressions for the wave functions are used and 
I > 0, due to the shape of the Hankel functions involved. In Refs. \2S\ I2TI] 
it was discussed that up to i = 2 this problem can be solved by calculating 
first the Fourier transform of the wave function in the k z direction and then 
getting the cross sections from the integral of the momentum distribution 
which is not divergent. On the other hand for any arbitrary £, total cross 
sections are well calculated using a numerical wave function while the check 
of peripherality can be done as discussed in the forthcoming Sec. 5. 

From the point of view of the impact parameter dependence, a similar 
comparison can be carried out. Fig. El displays the impact parameter profile 
of the same reaction. The 10 Be-n relative energy has now been integrated 
over, whereas the impact parameter is now along the x-axis. No cross-section 
is produced at low impact parameters due to the presence of the core survival 
probability \S ct (b)\ 2 . The maximum is reached at a distance near the sum of 
core and target radii, and then the cross-section decays with a long tail. The 
difference between the full wave function and its asymptotic approximation 
decreases for larger values of impact parameter, as these imply the neutron 
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Figure 4: Calculated impact parameter profile in nuclear elastic breakup of 
n Be on 208 Pb at 70 A.MeV. Full line: full-range bound state wave function. 
Dashed line: asymptotics. 



being far away from the 10 Be core, hence involving bigger values of r, for 
which the asymptotics are better. 

When nuclear breakup of the same system is looked upon, equivalent 
features come out. Fig. 0] shows a similar impact parameter profile for this 
11 Be reaction, but now the breakup is caused only by the nuclear interaction. 
Again, or perhaps more clearly in this case, asymptotic approximation and 
full wave function disagree the most for the lower half of the relevant impact 
parameter range, yielding almost identical results for the upper one. As 
expected, the decay of the profile with impact parameter is much faster now 
than in the Coulomb breakup case, what obviously accounts for the long- 
range character of the latter. In any case, these calculations show that no 
important effect is introduced by the approximation of the wave function by 
its tail. 



4.2 Validity of the dipole approximation 

As discussed in section |2J the phase given by eq. (fTHj) contains the Coulomb 
interaction to all orders in the multipole expansion. Thus the importance of 
high order multipole terms can be checked by comparing the results obtained 
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Figure 5: Calculations of differential cross-section for Coulomb breakup of 
11 Be on 208 Pb with respect to 10 Be-n relative energy in the centre of mass of 
the projectile. Full line: dipole approximation. Dotted line: full multipole 
expansion. 

with such expression with those provided by the first order term in its mul- 
tipole expansion, i.e. the dipole phase eq. (fT§|). In Fig. El such comparison 
is performed. As expected from earlier works, high order effects are of little 
importance in this case. In addition, they interfere destructively with the 
dominant dipole term, giving rise to a decrease in the single particle cross 
section, as has been also shown previously |2H1 EH] • 

4.3 All orders and nuclear effects 

As an extension of the work performed in [21], once the inclusion of all the 
terms in the multipole expansion of the potential has been achieved within 
this formalism, the effect of higher order terms in the strength of the po- 
tential is assessed next. In j2H] it was shown that a possible way to take 
into account this higher order strength terms is to use the sudden approxi- 
mation in which the (impact parameter divergent) first term is subtracted, 
and then a first order term, calculated in time-dependent perturbation the- 
ory, is added. This technique is now revisited with our full-multipole version 
expressions for Coulomb phase and its sudden approximation, eqs. ()18j) and 
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Figure 6: Differential cross section as a function of core-target relative energy for 
Coulomb breakup of 11 Be on 208 Pb at 70 A. MeV. Full line: all orders. Dotted line: 
first order time-dependent perturbation theory. 
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Figure 7: Differential cross section as a function of core-target relative energy for 
Coulomb breakup of 11 Be on 208 Pb at 70 A. MeV. Full line: all-multipole-all-order 
calculation for Coulomb breakup. Dashed line: including the nuclear interaction 
in the sudden approximation. Dotted line: nuclear breakup |28j . 
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Figure 8: Our calculation for the relative energy spectrum following 11 Be 
elastic breakup compared to full angular coverage measurements. Full line: 
theoretical calculation with spectroscopic factor 1. Dotted line: theoreti- 
cal calculation scaled with 0.6 to match the peak of the experimental data. 
Squares: data from j3] . 

()20|) respectively. Moreover, the full-range bound state wave function is used 
instead of its tail, although this effect has been above proved to be small. 
Fig. IH1 shows the result. Again, a first order approach seems to suffice to 
account for the vast majority of the physics involved. 

Finally, notwithstanding the fact the elastic breakup of n Be on Pb is 
well-known to be Coulomb dominated, the nuclear breakup are included in 
this all-order multipole-and-strength calculation, and the result showed in 
Fig. El Unsurprisingly, the inclusion of nuclear breakup makes very little 
alteration to the only-Coulomb result. 

4.4 Comparison to data 

Let us now turn our attention to the matter of how these calculations per- 
form when compared to experimental data. Figs. |H1 and El explore this issue. 
If the whole angular range measurements are taken (Fig. EJ), the spectro- 
scopic factor suggested to match the peak of the experimental data would be 
0.6. However, the shape of the distribution is poorly described, what shed 
some suspicions about the suitability of our method to calculate this partic- 
ular observable. On the other hand, Fig. E] shows the same comparison for 
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Figure 9: Forward-angle relative energy spectrum for elastic breakup of 11 Be. 
Full line: theoretical calculation with spectroscopic factor 1. Dotted line: 
theoretical calculation scaled with 0.5 to match the peak of the experimental 
data. Squares: data from jl] . 

forward- angle measurements. In order to match this requirement, the usual 
semiclassical relation between scattering angle and impact parameter has 
been used jl] , in such a way that a minimum impact parameter is given for a 
maximum scattering angle. We have taken this maximum impact parameter 
value as 20 fm. When this is carried out, the peak of the data is reached by 
the inclusion of a factor 0.5. Besides, the shape of our calculation agrees now 
much better with that of the entire experimental distribution. This is not 
surprising as both in the eikonal approximation and in perturbation theory 
we use a straight line trajectory which works at best in processes happening 
at high impact parameter, and where high-order, multi-step processes are of 
little relevance. A similar result has been obtained by the numerical solution 
of the neutron Schrodinger equation method of Baye and collaborators |2Z] ■ 
The analysis in this section has been performed for a typical halo nucleus 
at an energy at which experiments have been performed at various laborato- 
ries. Another interesting projectile to study [3H] would be 19 C whose initial 
state wave function contains also d-components and whose experimental neu- 
tron separation energy is not perfectly determined. Also, for this nucleus the 
numerical solution of the Schodinger equation is probably quite complicated. 
In the next section we will indeed study the sensitivity of our results to sepa- 
ration energies and to the initial angular momentum for heavy projectiles. As 
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far as the incident energy is concerned our conclusions are certainly valid at 
higher energies, while incident energies close to the Coulomb barrier cannot 
be treated by the eikonal approximation. 

5 Tightly bound systems 

For non-halo systems, this formalism can be used to investigate the gen- 
eral trend of nuclear and Coulomb breakup cross-section as a function of 
the separation energy. The projectile under study is now 34 Si whose valence 
d-orbital neutron has a separation energy S n =7.54 MeV. Experimental data 
already exist jSHj • From now on we shall include also stripping in our calcula- 
tions since this component is dominant for normally bound nucleons. Before 
showing the results for total breakup cross sections we wish to investigate 
the question of the peripherality of the reaction. 

As mentioned early on, nuclear breakup is calculated in the sudden ap- 
proximation, thus our formulae for elastic breakup and stripping cross sec- 
tions, after averaging over angular momentum projections, are 

a el = J db\S ct (b)\ 2 J dr\Mr)\ 2 \l-S nt (\b + r\)\ 2 (24) 

and 

<r« = J db\S ct {b)\ 2 j dr\Mr)? (l - \S nt (\b + r\)\ 2 ) . (25) 

Here S nt (\b + r\) — e lXnt ^ b+7 ^ defined by eq. In order to address 

peripherality, we study how different elements of the integral over r in above 
equations behave for a fixed value of impact parameter. We show them in 
Figs lTUl ITTl IT21 and El for core-target impact parameters of b = 12 fm and 
b = 15 fm, two values close to that where cross section is maximum. The 
terms |1 — S\ 2 and 1 — IS*) 2 have been projected onto the direction of r for 
simplicity. The asymptotic part of the wave function is also shown. This 
matches perfectly the wave function from 5 fm onwards. Similar figures for 
a much lighter system ( 12 Be+ 9 Be) where shown and discussed by Ref . |30] . 
It is clear from these figures that although the integral in eqs. (J24j) and (J25|) 
extends up to r=0, the integrand is concentrated at the surface. Furthermore 
smaller core-target separations are excluded by the experimental detection 
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Figure 10: Different factors entering the integration over r in the calculation 
of the elastic breakup cross-section, eq. ({21} for a core-target impact param- 
eter of b = 12 fm. Full line: radial wave function; dotted line: asymptotic 
form; dotted-dashed line: the factor |1 — S nt {b — r)| 2 ; dashed line: integrand 
of r-integration in eq. ()24|) magnified by a factor 50 to make it visible. 
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Figure 11: Same as Fig. El for an impact parameter of 15 fm. 
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Figure 12: Same as Fig. El but for stripping with a core target impact 
parameter of b = 12 fm. Dotted-dashed line is now 1 — \S nt (b — r)| 2 . The 
magnifying factor is now of 2000. 
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Figure 13: Same as Fig. [121 for an impact parameter of 15 fm. 
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condition of a surviving core. At the same time the target, four times heavier 
than the projectile, will suffer a negligible excitation and no recoil. 

Next, Fiefm shows the result for the breakup cross sections. In the 34 Si 
system, bound states with 2.5, 7.5, 10.5 and 14.5 MeV of separation energy 
have been artificially generated as s, d and f-waves, and subsequent breakup 
cross-sections have been calculated using a Pb target at 70 A. MeV. As Fig. 
El displays, for the three angular momenta chosen the nuclear breakup seems 
to reach a limit value of separation energy, where cross- sect ions depend just 
on nuclear sizes. On the other hand, Coulomb breakup continues decreasing 
as separation energy increases, to a point in which, even for a heavy target, 
becomes negligible. The curve suggests a dependence a CO ui oc exp (— ej/2). 

Finally on FigsEl and we show the core parallel momentum distri- 
butions after one neutron breakup (with S n = 8.07MeV, £=3 and £=2 re- 
spectively) of 46 Ar on Pb at 70 A. MeV. The eikonal total cross sections (full 
line) calculations show the typical shapes characteristic of nuclear breakup 
in the same way as previously seen for much lighter and less bound exotic 
projectiles on light targets. The corresponding cross sections are a st = 15.7 
mb and <7 e /=12.0 mb for the initial d-state, a st =lQ.O mb and cr e i= 11.9 mb for 
the initial f-state. The curves include elastic breakup and absorption which 
are also separately shown in the same figures by the dotted and dashed lines 
respectively. This last example shows that the deviations from the eikonal 
model seen in the experimental data of Ref . j3I] could be overcome by using 
a heavy target, for which the peripherality and no-target-recoil hypothesis 
would hold. For completeness we show also in Fig. El the experimental free 
n- 208 Pb cross sections and the calculated elastic, inelastic and total cross 
sections for the same system. We used the effective nucleon-nucleon interac- 
tion of Jeukenne, Lejeune and Mahaux (JLM) EZj, folded with the Pb 
density [33] to obtain the n-Pb optical potential. 

6 Conclusions 

In this paper we have presented a - mostly analytical - method to take into 
account the effect of the full Coulomb potential to all orders in the interaction 
in the semi classical scattering of two nuclei. The corresponding core- recoil 
effect for an exotic nucleus projectile described as a core plus neutron sys- 
tem, gives rise to the so called Coulomb breakup. We have then calculated 
Coulomb and nuclear breakup to all orders and checked the "peripherality" 
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Figure 14: Integrated breakup cross-section for a hypothetical Si beam 
against Pb at 70 A. MeV as a function of the neutron separation energy. Dif- 
ferent initial parameters: circles (squares) are for Coulomb (nuclear) breakup 
from an initial s-wave; diamonds (triangles) for Coulomb (nuclear) breakup 
from a d-wave; pluses (stars) for Coulomb (nuclear) breakup from an initial 
f-wave. Nuclear breakup is the sum of diffractive and stripping contributions 
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Figure 15: Model prediction for the core parallel momentum distribution of 
45 Ar cores from an initial £=3 state. The stripping (elastic breakup) part is 
the dashed (dotted) line. Total is the full line. 
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P„ (GeV/c) 

Figure 16: The same core parallel momentum distribution as in Fig. ED 
with an £=2 bound state. The stripping (elastic breakup) part is the dashed 
(dotted) line. Total is the full line. 
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Figure 17: Experimental total n- 208 Pb cross sections [42J and the calculated 
JLM jSZj elastic (dotted line), inelastic (full line) and total (dashed line) 
cross sections 
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of the reaction for light and heavy targets and for halo and normally bound 
neutron rich nuclei. The figures show that the effects of the asymptotic ap- 
proximation, higher multipoles and higher order in perturbation theory are 
accurate to within a few percent. 

Finally we have shown that, for the medium mass to heavy exotic nuclei 
that have started to be studied recently and will be further produced in the 
future, the Coulomb breakup effect on a heavy target becomes negligible with 
respect to the nuclear breakup part. This happens because of the "normal" 
separation energies ( >10MeV) and the large angular momentum of the va- 
lence orbitals i ^ 2. Therefore we have suggested to study these kind of 
nuclei by nuclear breakup reactions on heavy targets since all the conditions 
to apply the core-spectator model and the no recoil approximation for the 
target are satisfied. 
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